ON THE EIGENVALUE SPACING DISTRIBUTION FOR A 
POINT SCATTERER ON THE FLAT TORUS 
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Abstract. We study the level spacing distribution for the spectrum of 
a point scatterer on a flat torus. In the 2-dimensional case, we show that 
in the weak coupling regime the eigenvalue spacing distribution coincides 
with that of the spectrum of the Laplacian (ignoring multiplicties), by 
showing that the perturbed eigenvalues generically clump with the un- 
perturbed ones on the scale of the mean level spacing. We also study 
the three dimensional case, where the situation is very different. 



1. Introduction 

1.1. The Seba billiard. Point scatterers are toy models used to under- 
stand aspects of quantum systems for which the corresponding classical limit 
is intermediate between integrable and chaotic. In this paper we study the 
spectral statistics of point scatterers on the flat torus (a " Seba billiard" ) in 
the "weak coupling" regime. 

A point scatterer on the torus is formally given by a Hamiltonian 

(1.1) -A + a5 X0 , a CM. 

where A is the Laplacian, a denotes a coupling constant and xq denotes the 
position of the scatterer. Mathematically a point scatterer is realised as a 
self-adjoint extension of the Laplacian —A acting on functions which vanish 
near xq (see [3]). Such extensions are parameterized by a phase ip 6 (— it, tt], 
where <p = it corresponds to the standard Laplacian (a = in (jl.ip ). We 
denote the corresponding operator by —A XOjip . For (p = tt the eigenvalues are 
those of the standard Laplacian. For ip ^ it (a ^ 0) the resulting spectral 
problem still has the eigenvalues of the unperturbed problem, with multi- 
plicity decreased by one, as well as a new set of = {Aj} of eigenvalues 
interlaced between the sequence of unperturbed eigenvalues, each appearing 
with multiplicity one, and satisfying the spectral equation 



j=0 \ J 3 
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where cq = Yl'jLo ^xt+i anc ^ {^j} form an orthonormal basis of eigen- 
functions for the unperturbed problem: —Aipj = Xjipj- The eigenfunction 
corresponding to A G A v is the Green's function G\(x, xq) = (A + A) -1 ^. 
We denote the unperturbed eigenvalues without multiplicities by 

= {n = < m < • • • < rij < . . . } 

and call them the "norms" of the torus. Note that the perturbed eigen- 
values defined by (|1.2p are independent of the location xo of the scatterer, 
since in the case of the torus the sums Yl\j=n iV'j^o)! 2 = = n } are 

independent of xq. 

1.2. Spacing distributions. The perturbed eigenvalues {Aj } interlace with 
the norms {%} as follows 

(1.3) A£ < = n < A^ < m < • • • < \l < n k < ■ ■ ■ 

The nearest neighbour spacings for the norms and for the perturbed eigen- 
values are defined by 

(1.4) 6,: n r 5j := Aj +1 - Aj 
The mean spacing between the norms is defined by 

rij<x 

where 

(1.6) N{x) := #{j : n, < x} 

and likewise for the mean spacing between the new eigenvalues. 

Clearly 

(1.7) 

We define normalized nearest neighbour spacings by 
(1.8) 5 3 := * SJ : = -j-^— 

We want to determine the distribution of the normalized spacings 5j. 

Shigehara, Cheon et al. [9l I1Q|. [T2] identify two regimes in the semiclassical 
limit for a point scatterer in dimension 2: In the weak coupling regime the 
phase tp is fixed as A — > oo. In this regime the authors predict a Poissonian 
level spacing distribution for the perturbed spectrum. The strong coupling 
regime is when ip varies as A — > oo so as to satisfy: cotan^ ~ — ^logA. 
where they predict level repulsion. In most numerical studies of this prob- 
lem, it is the second regime that appears, due to a truncation procedure [8]. 
For an analytic study of this regime see [HJ [6] . 
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We deal with the spectrum of a point scatterer in the weak coupling regime 
(if fixed). We will show that the level spacing distributions of the norms and 
of the perturbed spectrum (if either exists) coincide. Since it is generally 
believed that the spacing distribution of the norms is Poissonian (if the torus 
is either rational, such as the standard torus M. 2 /Z, 2 , or generic irrational in 
a suitable sense [4]), that would imply that the perturbed spectrum is also 
Poissonian. 

1.3. Our results. We denote the differences between the old and new eigen- 
values by 

(1.9) dj := rij - A J > 

Since Sj — Sj = dj+i — dj, the normalized nearest neighbour spacings between 
the norms and the perturbed spectrum are related by 



(1.10) Sj -6J~- 

We define the mean difference of dj by 



dj+i - dj 



Xf<x 

We will show that the ratio between the mean difference dj and the mean 
spacing 5j vanishes: 

Theorem 1.1. For a point scatterer on a 2 -dimensional flat torus, 

(1.12) Trr^O, as x-^oo. 

As a consequence, since the differences dj > are non-negative, we deduce 
Corollary 1.2. Outside of a zero-density subsequence, 

(1.13) _L o, as j -> oo . 

That is the norms and the perturbed eigenvalues clump together generi- 
calljQ on the scale of the mean spacing. Therefore: 

Corollary 1.3. If the spacings Sj for the norms have a limiting distribu- 
tion, then so do the spacings Sj for the perturbed spectrum and the limiting 
distributions coincide. 

A similar result holds for hyperbolic surfaces if the point scatterer is 
placed in a generic position. We will not give the details here. 



^Recently Tudorovskiy, Kuhl and Stockmann [14] presented a heuristic argument that 
in the fixed regime the spacing distribution should be Poissonian by claiming the bound 
(|1.13[) holds individually, for all j. We are unable to verify this. 
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1.4. Dimension 3. The situation is very different for a three-dimensional 
torus T 3 . Let rjj be the perturbed eigenvalues of the point scatterer and rjj 
the unperturbed eigenvalues counted without multiplicity (the norms). The 
ordering is 

(1.14) r% < = r] < rjf < m < • • • < Vj < Vj 

As before we let dj := rjj — rf- and Sj = r/j+i — rjj. We denote by {dj) x 
the average of the spacings dj, and by (Sj) of the spacings of the norms, 
for r/J < x. 

Theorem 1.4. For the 3- dimensional flat torus, we have 

(1.15) ton #^ = -. 

Note that Theorem 11.41 does not give any information on the relation 
between level spacing distributions for the norms and for the perturbed 
spectrum. For an empirical study of the spectral statistics in dimension 3, 
see p]. 



2. Overview of the proof 

2.1. Our method. We derive Theorems 11.11 and 11.41 from the asymptotics 
as /3 — >• of the sum 

oo 
3=0 

To so, we approximate the sum by the difference 

oo 

of the heat traces of the operators — and —A, which we study in § [5] 
by using a trace formula which will be developed in §[3jH] in the 2-dim case. 
The 3-dimensional case is treated in § [6l 

2.2. A trace formula for the point scatterer on the torus. We work 
with a rectangular two-dimensional flat torus T 2 = M. 2 /2irCo, where Co = 
Z(l/a, 0) © Z(0, a) for some a > 0. Denote by C the dual lattice of Cq. The 
eigenvalues of the Laplacian on T 2 are the norms of the vectors of the dual 
lattice C (cf. section 2 in [7]). We denote the set of norms of the dual lattice 
vectors by 

(2.1) M = {0 < m < •••} 

and the multiplicity of an eigenvalue n G N is denoted by 

(2.2) r£ ( n ) = #{£ £ £ : |£| 2 = n}. 
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Recall that the perturbed eigenvalues {Aj} interlace with the norms {rij}. 
The ordering is 

(2.3) A£ < = n < Xf < n\ < ■ ■ ■ < Aj < rij. 

(That A£ < is given in U). 

We denote rij = p|, where pj > for j > 1, and A J = (pj) 2 , where pj > 
if j > 1 and 9/9q > (note that Aq < and p^ is pure imaginary). The 
spectral function 

1 °° f 1 1 

(2.4) S v (p) =-^ + E | — j " ^ | " co tan f 

has simple poles at the points p = ±pj and zeroes at the points p = ±pj . 
For a large enough and Qp = — a we will show that 

(2.5) S» = -i- log(ip) + i-£(p) + c(^) 

for c((/?) = ci — Co tan ^, where c\ is some real constant, and \D(p)\ <C CT 1. 

Let h be an even function which is analytic in a strip < a' for some 
a' > a and satisfies 

i%)i«(i + iK P ir 5 - 5 

for some 5 > uniformly in the same strip. We have the following general 
trace formula which we prove in sections El HI Let a > oo(<p) be sufficiently 
large. Then for all h as above, we have 

(2.6) 

"»<«»- 557^, P (logi P -2»c( V )) 

2.3. A Tauberian theorem. To prove Theorem 11.11 we will employ a 
Tauberian Theorem and reduce the problem to studying the asymptotics 
as \ of 

(2.7) !(/?) = J>e-^ 

i 

To study ^4(/3) we prove the following approximation (cf. (15. lj) and Lemma 

EH) 

(2.8) = -e-^} + 0(r 1/2 ). 
i=o P j=o 

We then use the trace formula (|2.6p with /i(p) = e _/3p2 to bound ^j^oi 6 "^^ ~~ 
e~^™ J } and obtain the following estimate which is the key result in the proof 
of Theorem ll.lt 
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Proposition 2.1. As (3 \ 0, 

(2.9) I ( £) = ^ e -/5AI«^_ 

2.4. Proof of Theorem 11.11 we will use Karamata's Tauberian Theorem 
(see e.g. [5]) which deals with the following situation: We say a positive 
function L(x) is slowly varying if L(kt) ~ L(i) as t — > oo for each fixed 
k > 0. We are given a non-decreasing function A{t) on M + such that the 
Laplace transform 

/>oo 

(2.10) A{p) := / e- pt dA(t) 

Jo 

converges for all (3 > 0. Suppose there exists two real numbers c > 0, u) > 
and a slowly varying function £(x) so that 

(2.11) A(f3) = {c + o(l)}p-"L(l/P), 0\O 
Then 

(2.12) A{x) = {c + o{l)}^^ y x^oo 

We apply Karamata's Tauberian theorem to the function 

(2.13) A{x) := dj 

\j<X 

which is non-decreasing since dj > 0. The Laplace transform A is 

(2.14) A{P) = Y,d 3 e~^ . 



Proposition 12.11 implies that A{j3) = o(l/(/3y/log -|)). Thus in Karamata's 
theorem, we may take oj = 1, L(t) = l/^logt, and c = to find 

(2.15) A(x) = o ( X ) , x-^oo. 



\/logx 

Therefore 

(2W) _ A{x) N{x) _ A{x) 

1 ' ' (Sj) x N(x) x x " V v 7 !^ 

as x — > oo, proving Theorem 11.1 



2.5. Three-dimensional tori. As in the 2-dimensional case, Theorem 11.41 
follows from the following proposition which we prove in section [6l 

Proposition 2.2. We have as (3 \ 

oo 1 

(2.17) Y. d ^~ H ' = TB + ° {rV4) - 

3=0 P 
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The key tools in the derivation are a trace formula (cf . Theorem 16. 3p and 
an approximation lemma (cf. Lemma 16 . 5|> . 

3. The trace formula 
We follow the same path as in |15j for a compact quotient T\H. 

3.1. Overview of the proof. Let T > be such that T ^ {pj} U {pj}, 
a > and consider the box 

B(a,T) = {p\\Zp\<a,\Mp\<T}. 

For a large enough and Qp = —a we will show that the spectral function 
(12.41) can be written as 



(3.1) S^p) = -i- log(ip) + ±-D(p) + c(<p) 

for c(ip) = c\ — Co tan ^, where c\ is some real constant, and \D(p)\ <C<7 1. 

Let h be an even function which is analytic in a strip |Qp| < a' for some 
a' > a and satisfies 

|%)|«(1 + |^|)- 5 - 5 
for some 5 > uniformly in the same strip. A contour integration gives 

(3.2) 2 £ %J)-2 £ h(p j ) = -±-[ h(p&p)dp 

We may rewrite (|3.2p as 

2fcO#)-2ft(0)+2 £ fc(pp-2 J]) 

0<pJ<T 0<p f <T 

(3 .3) i r iCT+T 



where 



vri 7_ i(T _ T p(log ip - 27rc(v9)) 



- / " T h'(p)log(l- . V{p) )dp + dB(T) 



i r / a r~\ \i -i<H-T 

8B(T) =— h(r)log 



(3.4) 



7T1 



S v (r) 



log(ir) — 2irc((p) 



-ia-T 



1 /' i,7 + T S" 
+ -/ Hp)^(p)dp. 

71-1 J-ia+T ^ 



Choose a sequence {T n } away from {pj} U {pj} such that lim n T n = oo. 
By use of the asymptotics (|3.1j) we show that the integral over the contour 
[—icr — T n ,—ia + T n ] converges absolutely as n — > oo. Since Weyl's law 
implies that both traces converge absolutely, it follows that lim„ dB(T n ) 
exists. The main step in the proof of the trace formula is to show that 
actually lim n dB(T n ) = for a suitable choice of a sequence {T n }. 
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In Lemma EL"T1 we construct a sequence {T n } which satisfies 

\S v (T n + iw)\ < e T*+ € . 

We then use this bound together with our knowledge of the existence of 
lim n dB(T n ), which holds in particular for a certain test function /15 with 
suitable symmetry properties (cf. Lemma [4.2p . We exploit the properties of 
this particular test function to bound log \ S^\ on average on the segments 
[— ia + T n ,T n ], namely (cf. Lemma FOl) 



(3.5) 



log 1 5^)1 dp 

i<r+T n 



which allows us to pass to the limit T n — > 00 and obtain the trace for- 
mula (ESI). 



Remark. We are unable to obtain an individual bound on log (S^p)) on the 
segments [—ia + T n ,T n ], but a bound on average suffices for our purposes. 

3.2. The Green's function on the torus. The free Green's function on 
R 2 is given by 

(3.6) g\(x,x ) = —K {ip\x - x \), A = p 2 , 

where Kq denotes the zeroth Bessel function. 
From the integral representation 

(3.7) K (r) = / exp(-r cosht)dt w=c ^ sht / % „ , 3?r > 

Jo J 1 yur — 1 

we obtain the following integral representation for the free Green's function 
on R 2 

I /-oo e -wp\x-x \ dw 

(3.8) gx{x,x ) = — , A = p, 9 ; p<0. 

2vr J 1 Vw 2 - 1 

We derive an integral representation for the Green's function on the torus 
T 2 by the method of images. Let Sjp < 0. We have 

k p (w; x, xo)dw 



1 f c 

(3.9) G x (x, x ) = ^2 9\(x + n,x ) = — / 

„ r- P J 1 



Vu^l 



where 

(3.10) k p ( W ;x,x ) = Y J e ~ iWPlX ~ X0+nl - 

Absolute convergence follows from the inequality (note uu > 1) 

(3.11) \k p (w;x,x )\ < e -«H*-*o+n| < 1+ ^ r £ (m)e^^™ 

The Bessel function has the asymptotics 

(3.12) K (z) = -log(z/2) -7 + o(l), z^O 
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where 7 denotes Euler's constant. Therefore the free Green's function has 
the asymptotics 

1 7 

(3.13) g x (x, xq) = -— logO|x - x \/2) - — + o(l) 

as x — )• xo- Thus we have the following asymptotics for the Green's function 
on the torus 

(3.14) G x (x, x ) = log(ip|x - xo|/2) - + C x + o(l) 
as x — > xq, where 

C\ = ^2 9x{xo + n,x ). 
nec\{0} 

3.3. In view of the spectral expansion of the Green's function G\ the spec- 
tral function may be written as 

(3.15) SJp) = lim {Gx(x,x ) - $tGi(x,x )} - c tan — 

x— >xq Z 

where 

\ - r £ (re) 

c = 1 + > 2 1 1 ' 

We may rewrite (13. 15ft as 

(3.16) S v (p) = ~ logip + i-P(p) + cfa) 
where 



(3.17) k(x)= e ~ ixH = E r ^ 

ne£\{0} O^meA/" 



k(pw)dw 
1 Vu> 2 — 1 



(3.18) 
and 

(3.19) c(<p) = --^P(-e i7r / 4 ) - c tan £ 

is a real constant. 

We have the expression 

(3.20) c(f) = c\ — Co tan — 
where 

cos( v / fw)e-v / f w ^ 



Vw 2 -1 



1 /" c 
(3.21) Cl = - Y, rc(m) / 

Lemma 3.1. For sufficiently large a > and 9p = —a 
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Proof. We have 



(3.23) | logiyo — 2ttc(<^)| > | log v^ 2 + i^p) 2 ~ 2ttc(^)| > log cr — 2tt|c(c^)| 
and 

(3.24) |P(p)| < / 1 )== < £ r £ (m) / = /(a) 

which implies for sufficiently large a > (in particular it is necessary that 
log a > 2ir\c((p)\) 

(3.25) „ Jg<gL_ < , ^> , „ < 1. 

| logip — 27TC(V?) I logfT — 2TT\C{tp)\ 

□ 

Let /i be an even function, analytic in a strip |3fp| < <ro, for some o"o > cr, 
which satisfies 

(3.26) |/ l (p)|«(l + |Kp|)- 5 -' 5 
uniformly in the same strip for some 5 > 0. 

Remark. We restrict ourselves to a smaller space of test functions here to 
simplify the presentation of our arguments. It is possible to obtain the trace 
formula for any test function with uniform decay \h(p)\ <C (1 + |5ftp|)~ 2_<5 . 

Let T > 0. Define the box 

B(a,T) = {p\ |»p| <T, \%p\<a}. 

Proposition 3.2. Denote by rij = p 2 , pj > 0, the eigenvalues without 
counting multiplicities. The new eigenvalues which lie strictly between the 
rij are denoted by Aj = (pj) 2 - We denote > Aq = (pg ) 2 where pg is purely 
imaginary and 3pQ > 0. Let a > Qpg and T > s.t. T ^ {pj}j U {pj}j- 
We have 

2h(p%)-2h(0) + 2 £ h( P J)-2 £ h( Pj ) 

0<pJ<T 0<Pj<T 

If . , ,s' 



(3-27) =-L / h(p)^(p)dp 

=- / +/ Mp)^(p)^ 

7" VJ-ia-T J-ia+T) \ 

Proof. By contour integration and symmetry. This is clear in view of the 
spectral expansion 



1 °° f 1 1 

(3.28) S v {p) = -j 2 + X>(n,-) - " c «tan 



□ 
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We may rewrite (I3.27|) as 

2h{pl) -2h(0) + 2 £ %p-2 ^ frfo) 
0<pJ<T 0<Pj<T 



(3.29) 



where 



(3.30) 



™ p(log ip - 27rc(y?)) 

-io-+T 



- / " log (l - . ) dp + dB(T) 



dB(T) =— 

7T1 



o / \ \ ~i — kr+T 

^(r)log' ^ J 



log(ir) — 2-7rc((p) 



-ia-T 



+ -/ Hp)^(p)dp. 



We have the following fact, analogous to Theorem 12 in |15| . 

Proposition 3.3. There exists an increasing sequence {ti} C M+ \ ({pj}j U 
{pjij) su °h that lim^oo ti = +oo and 

lim &B(fy) = 0. 

I— >oo 

Since the sums and integrals in (|3.29p (where we take T = tj) converge 
absolutely as ti — > oo, Proposition 13.31 which we will prove in section HI 
gives the general trace formula: 

Theorem 3.4. Let h be as (|3.26p . Let a > be large enough s.t. condition 
(|3.22p is satisfied. We have 

oo 

X>(P?) - K Pj )} 

3=0 

(3.3i) = ± r ia+o ° kp)*p 

27ri J-ia-oo PQ°g ip - 2ltc{ip)) 
J-ia-oa V 



log ip — 2irc(<p) 



dp. 



We call the first term on the RHS of (13.31R the "smooth term", and the 
second one the "diffractive term". 

4. Proof of Proposition 13.31 

We begin with the following lemma. 

Lemma 4.1. There exists an increasing sequence {T n } C \ ({pj}j U 
{pjlj) such that limn-i.oo T n = +00 and for —a < w < we have 

(4.1) |^(T n + hi,)| « £ T n 4+e . 
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Proof. We can choose an infinite increasing subsequence of Laplacian eigen- 
values {n fc(n) }„ such that n fc(n)+1 - n k{n) = p 2 k{n)+1 - p 2 k{n) > 1. This is be- 
cause the mean spacing between the norms {nj} is of size -y/log rij if the lat- 
tice C is rational and of size 1 if the lattice is irrational. Recall that between 
two consecutive eigenvalues n k ^ = ph n \ and n k (n)+i = pL n ) +1 there is ex- 
actly one new eigenvalue A£ (n)+1 = X 2 k(n)+1 and Xk(n)+i G (Pk(n), Pk(n)+i) C 
1R+ is a zero of the function S v (p), whereas p k ( n ), Pk(n)+i are singularities 
of the same function. 

So we may choose an infinite sequence 



( i 



(4.2) T n 



2(Pk(n) + Xk(n)+l), if |Xfe(n)+l _ Pk(n)\ > \Xk(n) ~ Pk(n)+l\ 

^(p k (n)+i + Xk(n)+i), otherwise. 



with \p k ( n) - Pk(n)+i\ > \pk(n) + Pfc(n)+il ~ T n 1 . Note in particular that for 
all pj € R+, 



(4.3) 



\pj 



T n \ > -g\Pk{n) - Pk(n)+\ 



»T T 7 1 . 



Let p n {w) = {T n + iw) 2 , w € [— a, 0]. We have for any e > 



(4.4) 



\S*(T n + iw)\ «J> £ (n j 

j=0 



n ; 



Pn(w) 



<C e |i-^vH|J^ 



II; 



o l n i - MjvHIK 



where we have used the bound rc{n) <C e n e . Fix a £ (e>l)- We split 
the sum into a central part satisfying m( w( z^_ a ^ \rij — pn(w)\ < n" and a 
corresponding tail. For convenience we let I n {nj) = inf^^g] \ n j — Pn(w)\- 
The first sum is estimated by 



(4.5) 7 ™K')<S 



a \rij -p n (w)\\ nj -i| 



<#{i I < } max 



sup 



Now if rij > T% then I n {jij 



(4.6) 

Let 
(4.7) 



'«(» ; )<«° to6[-<r,0] I l n i - Pn(w)\\nj - i| 

rij - T%. It follows 

#{j | I n (nj) < n°} 
<#{j | < T n 2 } + #{j | n 3 - nf < T 2 }. 

C(a) = #{j\n J <2 1 /^} 
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and observe that nj > 2 1 ^ 1 ~ a * ) implies nj" 1 < i. So nj > 2 1 /( 1 - Q ) together 
with re,(l — n" _1 ) < T 2 implies 



(4.8) < 2nj-(l - n^ 1 ) < 2T 2 . 

Hence 



#{j|n,(l-nJ- 1 )<T n 2 } 
< #{j | < 2 1 /d-«), n ,(l - n"" 1 ) < T 2 } 

(4.9) + #{j | nj > 2 1 /d-«), % (1 - n"- 1 ) < T 2 } 
<C7(a) + #{i|2 1 /(i-«)<n J <2r 2 } 

« Tl 

It follows that 

(4.10) #{j | I^n,) < nf} « T 2 . 

By the same observations as above we see that I{rij) < n°- implies nj < 
max{2 1 /( 1 -°),2T 2 }. Also for any j > we have (cf. flOD ) 

|Pj - r n | > - Pk(n)+l\ > ^nT 1 

which implies 

\rij - Hn{w)\ = \Pj ~ (T n + i^) 2 | = \Pj ~T n - iw\\pj + T n + \w\ 

(4.11) >\ Pj -T n \(p j + T n ) 

> 1. 

Since \rij — i| > 1 we have 

(4.12) max sup {-. p-rr. -1 < T*. 

/nK)<n? ^ e [_ CTi0 ] I Pi - Mn(™)IK - i| J 

The tail can be bounded as follows 

™J < n j 

^ \rij — u n (w)\ In,- — i| ~~ ^-f In,- — i| 

I n (n i )>n a 1 3 rnK J 1 I n {nA>n a 1 J 

(4.13) ' 3 ' ' 



< 



oo n e-a 



L * \Ti ■ — 1 



\n,j — i 

3=0 1 3 



Finally note that \fj, n (w) - i| < T 2 . □ 
Recall 



dB(T n ) =- 

(4.14) 



V(r) 



h{r)l ° g V 1 log(ir) - 2™^) 

+ -/ h{p)-^{p)dp. 
vri J- ia+Tn £>tp 



-iv+T n 



-\o-T n 
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We know that linin^oo B(T n ) exists (for any test function h which satisfies 
the uniform bound \h(p)\ <C (1 + |Kp|) -2 ~ 5 in the strip \^p\ < a - the 
decay which is required by Weyl's law to ensure that the trace converges 
absolutely) and we want to prove that the limit is zero for any test function 
which satisfies the uniform bound 

(4.i5) i%)i«(i + |K P ir 5 - 5 

in the strip |9fp| < a. 

For the first term we have, in view of < /(c) along Qp = —a, 



(4.16) 



log 1 



log(a ± iT n ) - 2ttc(<p) 

which implies that this term vanishes as n — > oo. 
For the integral an integration by parts gives 

ria+T n gl 

h(p)^(p)dp = [h(p)logS tp (p) 



1 



logT n 



(4.17) 



i<T+T n 



icr+T n 



h'(p) log S v (p) dp. 



To see that the first term vanishes as n — > oo, observe that the identity 
ATOP and the bound (ET21j) imply 

|log,S(±i(7 + r n )| =|log5(-ia=FT n )| 

(4.18) < |log|S(±i<r + r n )|| + |argS(±k7 + T„)| 

= io g io g r n + o(i), 

where we used | arg S(±ia + T n )\ <C 1 as n — > oo. Similarly we see 



io-+T n 

h'(p) log S v {p)dp 

-i<r+T„ 

We have the calculation 

ri<x+T„ 



icr+T n 



io-+T n 



h'(p) log \SJp)\dp + 0{T~ 



h'(p) log\S v (p)\dp 



(4.19) 



p-*~p 



—icr—Tn 



h'(-p)]og\S v (-p)\dp 



T„ 

-ier+T, 



h'{p) log \S v {p)\dp 
h'{p)log\S v (p)\dp 



where we used S v (p) = S (fi (p), and so the term 



/ia+T n 
h'{p)\og\S v {p)\dp 
ia-+T n 

converges to a limit as n — > oo. 



iv+T n 



{h'(p) + h'(p)}log\S^p)\dp 
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To obtain the result we require two lemmas. The first lemma constructs 
an even test function which is analytic in a strip and the real part of whose 
derivative satisfies a certain polynomial lower bound in T n on the line seg- 
ment [— icr + T n ,T n }. 

Lemma 4.2. Choose o~q > a . Let 

-1 



(p 2 + ° 2 o) 2 ' 

We have for t G [— ctq, 0] and for all sufficiently large n 



(4.21) 

as n — ^ oo. 
Proof. We have 



3tti 5 (T n + it) = \Vtti 5 (T n + it)\ » 



7^5 
- 1 n 



h' 5 (p) 



4p 



(p 2 + * 2 o) 3 ' 

Let t £ [— oo,0]. A simple calculation gives 

4(T n +it)(T*-t 2 + al-2iT n tf 



\Stti 5 {T n + it)\ 



(4.22) 



--ft 



4(T n + it)(r n 2 -t 2 + qg-2iT n t) 3 

((2* - tf» + og)* + 42**)* 



- 1 - II 



as n 



oo. 



□ 



The second lemma gives a bound on log | | averaged along the line seg- 
ment [-icr + T n ,T n ]. 

Lemma 4.3. We have the following bound 



(4.23) 



-io-+T n 



log | Sy,(p)| dp 



<C T 5 



Proof. We know there exists a constant c > such that for all n and w 6 
[—0", 0] we have 

|S^(T n + ry;)| <cT n 5 . 

In Lemma f4. 2 1 we prove the existence of a test function /15 which is analytic in 
the strip \^sp\ < a, satisfies the uniform bound |^s(/5)| <C (1 + l^pl) -4 in this 
strip and in addition h^(p) = h^{p) and fth' 5 (T n + iw) = \fth' 5 (T n + iw)\ » 
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T"? v We thus have 

nil) 



log\SJp)\dp 



\<T+T n 



(4.24) 



< T- 5 f - log(c- 1 r- 5 |5 v (T n + \w)\)dw + 0(T~ 5 logT n ) 

J —a 

r° 

« - / m' 5 {T n + iw) \og{c- l T-^\S^{T n + iw)\)dw 

J —a 



+ 0(T~ 5 log T n ) 

«e 1 

because \h§{p)\ <C (1 + |3f?p|) -4 uniformly in < a and therefore 



lim 

n— >oo 



icr+T n 



m' 5 ( P ) \og \s v (p)\d P 



exists. 

We obtain 



□ 



i(T+Tn 



< T, 



-5 



h'{p) ]og\S v {p)\dp 

l-ia+T n 

in view of the identity (|4.2Q|) . We also used that by Cauchy's theorem the 
analyticity and decay of h in |9fp| < o"o, where o"o > a, imply the analyticity 
of h' in < cr and the uniform decay 

|/ i '(p)|«(l + |^p|)- 5 - 5 

in the same strip. It follows that 

lim dB(T n ) = 

n->oo 

which proves Proposition 13.31 



5. Proof of Proposition 12.11 

5.1. We want to apply the trace formula in order to obtain information 
about the average spacing between new eigenvalues and old eigenvalues. 
Let h(p) = e~@ p , for small j3 > 0. Upon dividing through by /3 we can 
rewrite the l.h.s. of the trace formula (I3.3ip as 



(5.1) 



3=0 



-/3A? 



3=0 

oo 

: ^^e-^ + 0(/3- 1/2 ) 

3=0 



where dj = nj — Aj > 0. The last line follows from the following lemma. 
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Lemma 5.1. We have the bound 

> 1 _ e -Pdj 



(5.2) 



j=0 



e-W I 1 



« r 1/2 



Proof. For x > we have the inequality 



l-e~ x 
< 1 < x 



and the bound dj <C nV 4 for j > 1 (cf. the greedy algorithm in [7], p. 7). 
It follows 



£4 

3=0 



e~^ I 1 



1 _ e -Pdj 



J'=0 



2 -£A V 

e j 



(5.3) 



j=0 



n V 2 ie -^ +/3e -/3A- 



and the bound N !f (x) <C x permits us to bound the sum by the following 
integral: 



(5.4) 



°° POO 

E n ]+i e ~^ « Z 3 / ^ 1/2 e-^dx « /T 1/2 - 
i=o ^° 



□ 



5.2. The smooth term. We have the following bound on the smooth term. 
Proposition 5.2. As f3 \ 



(5.5) 



1 

2^ 



— icr+oo 



-/V 



rip 



p(log(ip) - 2vrc(v?)) 



1 



Proof. Denote by Cg the contour following a semicircle in the lower halfplane 
centered at the origin of radius 5, where e 2nc ^ > 5 > 0, starting from —S 
and finishing at 5. By shifting the contour across the pole at p = —ie 2nc ^ 
to the real line we obtain 



(5.6) 



1 
2vri 



— lcr+OO 



+ 



2vri 



+ 

C 4 



dp 



J /9log(i/9e- 2 ^))' 



Note that the integral over the semicircle vanishes as 8 — )■ 0. 
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We may pick the branch of the complex logarithm in such a way that 
arg(x) = 7r/2 if x < and arg(x) = 3ir/2 if x > 0. Then for real p ^ 

1 1 

log(ipe- 27rc ^)) ~~ log(|/9|e- 27rc (^)) +i(f + arg(/)e- 27rc (^))) 
(5 ' 7) _log(|p|e- 2 ^)) - i(f + arg(pe- 2 ^))) 
log 2 (|yo|e- 27rc ^)) + 7T 2 /4 

and it follows that 

1 f e-Pp 2 dp 

= _J_ f e -/V f +arg( P e- 2 ^)) 

2vr 7 R \(_ 5 ,5) / o(log 2 (|p|e- J M*0) + 

1 f 00 e-^ 2 dp 



2 



(5.8) 



/<5 p(log 2 (pe- 2 ^)) + 



Since 



2 J e -27vcM S r (l g 2 r + z j-) 

2 y-27rc( V )+log<5 i 2 + X 

1 r 00 e~^ cMe2t dt 

- / 5 as o — > 0. 

2 *2 + ni 



-oo t 2 + ^ ~ 2 ' 



we obtain in view of (|5.6|) 

1 e -Pp 2 dp 



^J-ia-oo plog(ipe- 27rc ^)) 
( } = pe ^ M _ } 1 r» l-e-^ 4 ^ )e2t dt 

+ 2 7-oo i 2 + ^ 

Let 7 = e 47rc ^/3. We proceed by dividing the integral on the r.h.s. into 
two integrals over the ranges (— oo, log7|) and [^^| log7|,oo) for some 



small e > 0. 

We then bound the first integral as follows 

(5.10) / ^ « 7 / e 2i dt = i 7 e 

J -co t H — j- J —oo 

where we note that 

|1- 6-7^1 <7e » 
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because t < ^rpl log 7) implies je 2t < Y- 
We bound the second integral by 



(5.11) 



1 



|l-e-T e \dt roc 

< 



t~ 2 dt 



2fe|i og7 | t 2 + ^ fe|iog 7 | " (l-e)|log 7 |' 



□ 



5.3. The diffractive term. We continue with the bound on the diffractive 
term. 

Proposition 5.3. Let h(p) = e~P p2 . As f3 \ 



(5.12) 



/>'(/» loft (1-^4),//, 
logip 



log^ 



Proof. (I3.22p allows us to estimate 



(5.13) 



^(-icr + a) log ( 1 - , P( ia + S/ „ ! ,1s 



00 
00 



logi(— kr + s) — 2irc(ip) 



< / |/i'(-icr + s) 



\V(-ia + s) 



-<is. 



logi(— id + s) — 27rc(ip)\ 
We have 

(5.14) \h'(-ia + s)| = 2/3\ - ia + s| | e"^^^^ 2 1 < 2/3(ct + IsDe^ 2 "^ 
and 

r oo e -a^w dw 



(5.15) |2?(-iff + s )| < ^ r c \ 



rn) 



l y/w 2 — 1 



/(*) 



Iog(ff + ia)| > ±log(cr 2 + s 2 ). 



and finally 
(5.16) 

We continue our estimate as follows (recall a > max{l, e 2wc ^}) 

/ \h(-ia + s)\- — — 
/ „ I ioe i( - 

(5.17) 

<8{3e^ 2 f(a) 



-ds 



logi(— id + s) — 2irc((p)\ 

(a + s)e-^ 2 ds 
i log(u 2 + s 2 ) - A-Kc(ip) 

and the integral is bounded by 

r°° (a + s)e~P s2 ds a 

L log(a 2 + s 2 ) - 4irc(ip) < (21og(j-47rc(^))/3 1 /2 J Q 



(5.18) 



e~ w dw 



+ 



se 



-? s2 ds 



log((J 2 + S 2 ) — 4:7Tc(tp) ' 
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Let 7 = e 47rc( ^/3 and £ = e~ 27rc ^V. We bound the second integral on the 
r.h.s. as follows 

r°° se-^ds _ 1 f 00 we- w2 ds 

Jo log(<r 2 + s 2 ) - 4ttc(^) ~/3 7 log(cj 2 + ^) - 47rc(v?) 

27 



1 P 2 

ttt / we ™ ciw 

>g£ 2 Jo 



oo 

2 



(5.19) ^° log(£ 2 + ^) 

1 

+ — — r / we~ w "dw 

for small e > 0. The last line follows, since for w > 7 e we have (assuming 
7 < 1) 

logger 2 + w 2 ) 21og w + log(l + ja 2 /w 2 ) 



(5.20) 



Since 



we have 



logi L " log(i) 



> 1 - 2e + O 



/ 7 we"™ 2 ^ = 0( 7 2e ) 
■/ o 



yi |logi(-i(7 + S )| log(i) 



□ 



6. The three-dimensional case 

The three-dimensional case is very different. Consider the three-dimensional 
flat torus T 3 = M 3 /27r£ 3 , wh ere £g = Z(l/a6,0,0) 0Z(O,o,O) ©Z(0,0,6) 
for some a, 6 > 0. Denote by C 3 the dual lattice of Cq. The eigenvalues of 
the Laplacian on T 3 are the norms of the vectors of the dual lattice £ 3 . We 
denote the set of norms of the dual lattice vectors by A^ 3 and the multiplicity 
of an eigenvalue n £ M 3 is denoted by 

(6-1) r3 (n) = #{£G£ 3 :|£| 2 =n}. 

Let tjJ be the perturbed eigenvalues of the point scatterer on T 3 and rjj 
the unperturbed eigenvalues counted without multiplicity (the norms). The 
ordering is 

(6.2) r% < = m < r? x < 771 < • • • < V J < nj. 

Our main result is the following. Let dj := r)j — rf- > 0. 
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Proposition 6.1. We have as f3 \ 



oo 



(6.3) E^ e ^ r = ^ + °(^ 3/4 )- 
3=0 P 

6.1. The Green's function in three dimensions. The free Green's func- 
tion on R 3 is given by (cf. [16], p. 842, eq. (3.4)) 

e -ip\x~x \ 

(6.4) g ri (x,x ) = — r, 3?p< 0, p 2 =r]. 

Air\x — xo\ 

We periodise to obtain the Green's function on T 3 . 

I e -ip\x-x +n\ 

(6.5) G„(x,x ) = — > 1 r 

v 7 A y 4vr ^ \x-x + n\ 

n£C 3 



In particular the deficiency elements are given by 



(6.6) G±i{x,x ) = ^~Y1 



\X — XQ + 7l\ 

COS 



\x-x +n\ 
V2 



4tt ^— ' \x — xq + n 

ne£ 3 



• > 2 

where we note that ±i — ' r 1 ~ 1 



The spectral function is given by 



where 



.7) S$(p) = lim (G x - 3tGi)(x,x ) - tan| = — + D$(p) 



(6.8) ^(p) = -tan| + i - 7 = + -E, 3 (n) 



r3(n) 



V i i ^ . , e p - e V2co Hv% 



J? 



where the first and second terms in (|6.7f) and (I6.8P come from the regulari- 
sation (where h = x — Xq) 



e 

lim 



h^O Aw\h\ 

(6.9) 1 - W + 0(H 2 ) - (1 - §)(! + 0(H 2 )) 
= lim rm 



4vr 4tt\/2' 

6.2. The trace formula in three dimensions. We require the following 
lemma. 

Lemma 6.2. For sufficiently large a > \^Pq\ we have for Qp = —a 
6.10 1 , Y n < 1. 

\p\ 
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Proof. For Qp = —a it can easily be seen from (j6.8[) that \D^(p)\ = O(l). 
Furthermore | yo | > a. So (I6.10P certainly holds for sufficiently large a. □ 

We have the following trace formula for a point scatterer in three dimen- 
sions. 

Theorem 6.3. Let h be as above and a > \^sp"q\ large enough such that 
(|6,10p is satisfied. Then we have 

oo 

£>(pp - h( P ,)} 

2m J-ia-oc V P J 

Proof. Following the argument in the proof of the trace formula for two 
dimensions we obtain for a > |9po I the trace formula (an analogue of Krein's 
famous trace formula) 

00 1 r— icr+oo 

(6.12) J>(PJ) " K Pj )} = — / h!{p) log S 3 (p)dp. 

j=0 J-ia-oo 

In view of (|6.7p we rewrite the r. h. s. as 



1 

2^i 



ti(p) log S 3 (p)dp 

icr— OO 
icr+oo 



us) -si/^MM^]* 



+ 7T1 ^(P)log 1 + —)dp 



27ri J-ia-00 " V p 

and the first term can be evaluated by integration by parts and shifting the 
contour 

-icr+oo 



1 r~™+°° h( P )dp 



.14) 



27ri V-i^oo p 
1 f h{p)dp 1 /" h(p)dp 



ZniJcs P 27ri 7m\(-5,5) P 



=3^(0) 

where we recall that for some small 5 the contour C5 denotes the lower 
semicircle connecting —5 and 8 on the real line. □ 
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6.3. Proof of Proposition 16.11 Let h{p) = e ^ p . In this case the trace 
formula gives us 



.15) 



oo 

£ 

3=0 



I 



-fa* 



-fa 



} 



-lcr+oo 



pe 



-fa 



log 1 + 



dp 



and in view of (|6.1U|) we have for Qp = —a the bound 



(6.16) 

which implies 

(6.17) 

and therefore 
(6.18) 



log 1 + 



AmD$(p)\ „\B%{p)\ 



\P\ 



-ICT+OO 



la— oo 
— ic+oo 



-fa 



log 1 + 



p 



\dp\ 



-^\\m P )\\dp\ 



dt = O 



^ oo 



e H 



-fa 



} 



j=Q 



+ o 



An analogue of the greedy algorithm is required for the proof of Proposi- 
tion We state this clS Si lemma. 



Lemma 6.4. Let dj := rjj 
(6.19) 



We have the bound 

1/8 



dj < i)j 

Proof. Recall that each rjj is of the form q(m, n, k) = am 2 + bn 2 + ck 2 for 
real numbers a, b, c > and integers m, n, k. We need to show that for each 
j we can pick m, n, k such that 

\r]J - q(m,n,k)\ < r/j^ 8 

where the implied constant does not depend on our choice. Let si = rf- — 

am 2 . Let m = ^JrjJ/a and observe that si = if- — am 2 <C {jfj) 1 ^ 2 < V^ 2 ■ 

Now let S2 = s\ — bn 2 and choose n = \J s\/b so that S2 = s\ — bn 2 -C s\^ 2 . 

\fs2~fc such that S2 — ck 2 <C s^ 2 . With the above 



Finally choose k ■ 
choices of m,n,k we have 

\r]f - q(m,n,k)\ 



S2 



ck 2 « 4 /2 « Sl /4 « J' 



□ 
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The following lemma implies Proposition 16,11 
Lemma 6.5. We have the following identity 



oo 



(6.20) d i e ~^ = \ ~ e ~^ J > + °(/ 3_3/4 ) 

Proof. It is sufficient to prove the bound 

(6,!) g^T( 1 _i^)« r .v 1 

For x > we have the inequality 

1 - e~ x 
< 1 < x. 

x 



It follows from the inequality and Lemma [67 

j=o v y j=Q 

oo 

(6.22) «/? ^ v] /4 e-^ + f3e~^ 

3=1 
oo 

3=0 

and the bound N v {x) -C x 3//2 permits us to estimate the sum by the following 
integral: 

°° />oo 

(6.23) [3 ^ ^fie-^ « /? / x^e'^x^dx « /T 3/4 - 

■ n " 



j=0 



□ 
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